J/ip cc production in e + e and hadronic interactions 



A. B. KAIDALOV 

Institute of Theoretical and Experimental Physics, B. Cheremushkinskaya 25, 
Moscow 117259, RUSSIA 
E-mail: kaidalov@heron.itep.ru 

Predictions of the nonperturbative Quark Gluon Strings model, based on the 1/N- 
expansion in QCD and string picture of interactions for production of states con- 
taining heavy quarks are considered. Relations between fragmentation functions 
for different states are used to predict the fragmentation function of c- quark to 
J/i/i-mesons. The resulting cross section for J/-0-production in e+e — -annihilation 
is in a good agreement with recent Belle result. It is argued that associated pro- 
duction of cc states with open charm should give a substantial contribution to 
production of these states in hadronic interactions at very high energies. 



Investigation of heavy quarkonia production at high energies provides an 
important information on QCD dynamics in an interesting region of intermedi- 
ate distances from l/m,Q to tqq, where tuq is the heavy quark mass and tqq 
is the radius of a heavy quarkonia state. For c and b-quarks this is the region 
0.05/m < r < lfm. In this region both perturbative and nonperturbative 
effects can be important. Production of J/^-mesons is studied experimentally 
in e + e _ -annihilation, jp, hp, hA and AA-collisions. Analysis of hadronic in- 
teractions show that the simplest perturbative approach (color singlet model) 
1 does not reproduce experimental data 2 . This observation lead to an intro- 
duction of the color octet mechanism 3 of heavy quarkonia production. In this 
approach a set of nonperturbative matrix elements is introduced, which is de- 
termined from a fit to data. A characteristic prediction of this approach is a 
large transverse polarisation of J/ip and ip' at large transverse momenta 4 is 
not supported by the Tevatron data 5 . 

A new mystery to the problem of heavy quarkonia production has added 
recent result of Belle Collaboration 6 on a large production of J/^-mesons with 
charmed hadrons. The observed cross section at y/s — 10.6 GeV is an order of 
magnitude larger than theoretical predictions 7 , based on perturbative QCD. It 
is interesting that at this energy an associated production of J /ip with cc-pair 
is the dominant mechanism of J/ip production 6 . 

In this paper a nonperturbative approach, based on 1 /N-expansion in QCD 
and string picture of particle production is used for a description of heavy 
quarkonia production at high energies. The model based on this approach 
(the Quark Gluon Strings model (QGSM) 9 ) has been successfully applied to 
production of different hadrons at high energies. It has been also used for 
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Figure 1: Diagrams for J/ip production in e+e annihilation. 

description of inclusive spectra of hadrons containing heavy (c,b) and light 
quarks n^ 12 * 13 . In QGSM the fragmentation functions, which describe tran- 
sitions of strings to hadrons in many cases can be predicted theoretically 9 ' 10 
and are expressed in terms of intercepts of corresponding Regge trajectories. 
We will show that the model naturally leads to the cross section of J/tp pro- 
duction in e + e~ annihilation consistent with the Belle result. Estimate of the 
contribution of the same mechanism in hadronic interactions indicates that it 
can be important at energies y/s > 10 2 GeV. 

Let us first discuss heavy quarkonia production in e + e~ collisions. In these 
reactions cc-pair is produced directly by a virtual photon. However a probabil- 
ity of transition of such a state at high energies (far above threshold of charm 
production) to J/tp is very small. A simplest diagram of QCD perturbation 
theory (Fig. la) corresponds to a transition to a white cc state with relative 
momentum characteristic to J/tp by emission of two hard gluons. This cross 
section is suppressed at high energies by a factor 4m 2 /s and at yfs = 10 GeV 
constitutes 10~ 3 of the total cc cross section 7 . 

J/tp production in association with extra charmed pair (Fig. lb) does not 
have this suppression, but contains a smallness due to production of this pair 
and a high threshold of the processes. At high energies this mechanism can 
be considered as a fragmentation of c(c) to J jip. Calculation in the lowest 
order of QCD perturbation theory 7 shows that this mechanism is important 
at energies yfs > 50 GeV, but at yfs — 10 GeV is smaller than the mechanism 
of Fig. la by an order of magnitude and is about 0.07 pb. This is in a sharp 
contradiction with Belle result: a(J/ip cc) = 0.87t°; 2 g ± 0.17 pb. 

Note that for states of comparatively large radius like J / ip and especially ip' 
or Xc a nonperturbative fragmentation can be important. Thus I shall estimate 
a fragmentation of c(c) into J ftp using the nonperturbative model mentioned 
above. In this model particle production is described in terms of production 
and fragmentation of quark-gluon strings. A behaviour of the fragmentation 



2 



functions is determined in the limit z — * 1 from the corresponding Regge limit 
and is expressed in terms of Regge intercepts cti(0) 9 ' 10 . The fragmentation 
function of c-quark to J/ 1\) in this model is written in the form 10 

Df =a^z- a ^°\l- z )- a ^ +x (1) 

where a^( ) is an intercept of the J /ip Regge trajectory, which is known from 
analysis of data on spectrum of cc states and a nalysis of inclusive spectra 
of charmed particles (see below), A = 2a' D ,p 2 LD « 1. Thus this fragmenta- 
tion function is characterized by one constant a^. In order to determine this 
constant we will use a relation between fragmentation function of c-quark to 
J '/ip and fragmentation function of a light quark to Z?(D*)-meson in the limit 
z — > 1. According to rules formulated in refs. 10 ' 15 both functions have the same 
behavior on z: (1 — ^(""^(OH*) as z — > 1 and differ only by a kinematic factor 
related to mass difference between J/ip and Z?(D*)-meson 

-« D .(o)) 

(2) 



The quantities soi will be determined below. 

Now we shall find the fragmentation function in the limit z — * 1. In 
this limit it is related to the fragmentation function of a light quark to 7r meson 
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(1 - z ) 2 ( Q p(°)- a B.(o)) (3) 

where a p ,ctD*(0) are intercepts of p and D* Regge trajectories. They are 
related to a$ (0) by the following equation 14 

a p (0) + a^(0) = 2a D .(0) (4) 

I shall use the following values for these intercepts: a p (0) = 0.5,0^(0) = —2 
and olb* (0) = —0.75 in accord with eq.(4). An uncertainty in the value of 
a^,(0) discussed in rcf. 11 is eliminated at present by experimental data on 
inclusive spectra of charmed hadrons in hadronic collisions. 

The gamma functions in eq.(3) appear from Regge residues of the corre- 
sponding trajectories, which were chosen in accord with dual models are in 
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a good agreement with data on widths of hadronic resonances 16 . The cou- 
pling is assumed to be universal (with an account of SU(4) and heavy quark 
symmetry) . 

The quantities Sot entering in eq.(3) can be easily calculated using formulas 
and parameters of ref. 15 

( s «c)2a D ,(0) = ( s u«)a p (0)( s Z5fl^(0). 

4 U = *™l± = 1 GeV 2 ; s° D = (m c± + m ul _f 

With m uA _ = 0.5 GeV and m cA _ = 1.6 GeV 15 we obtain Q g c ) = 3.57 GeV. 
Using these values for Soi m eq.(3) and m 2 ± = 0.18 GeV 2 , m 2 D± = 5 GeV 2 , 
the fragmentation function £>J + = 0.44 9 we obtain the function D® + at z — > 1 
in the form 0.01(1 — z)^ Q *' ' +A '. This value is in a reasonable agreement with 
phcnomcnological studies of charmed particle production in hadronic interac- 
tions in the framework of QGSM 12 ' 13 . 

The value of Sq^, in eq.(2) can be calculated in the same way with the 
substitution s^ 5 — » Sq D = 6.72 GeV 2 . Finally we obtain from eq.(2) 

Df = 0.05 (1 - z)(- a ^°) +A >; z -> 1 (7) 

Thus a^, = 0.05. 

At asymptotic energies s — > oo cross section for J/ip production in e+e~ 
annihilation is equal to 

a^=2a cB [ Df(z)dz (8) 
Jo 

factor 2 in eq.(8) takes into account J/ip production by both c and c quarks. 
At energy i/s ~ 10 GeV there is an extra suppression due to phase space 
corrections for production of a heavy state. We estimate it by introducing an 
extra factor 7 = yl — 4M§ /M 2 B to eq.(8). Distribution in M 2 B is related to 
the z distribution. It has a maximum at M 2 « 0.27s. For energy of Belle 
experiment the correction factor 7 = 0.7. Thus we obtain the following cross 
section for J/ip cc production at ^/s = 10.6 GeV — 1.2 pb. This value is 
in a good agreement with Belle result 6 and is much larger than perturbative 
QCD prediction 7 . An estimated uncertainty in the value of cross section due 
to possible variation of quantities SQi,rrii± and cti(0) is about 50%. 

Let us consider now J/ip production in hadronic interactions. In the ap- 
proach based on l/7V-expansion 8 the main diagrams for particle production 
correspond to two-chain configurations, shown for pp-interactions in Fig. 2a 
9 . They can be considered as production and fragmentation of two q — qq 
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Figure 2: Diagrams for J/ip production in pp- interactions. 

strings. It is important to emphasize that production of one cc-pair together 
with light quark pairs in this approach always leads to an open charm produc- 
tion (Fig. 2b) and J/ip in this case is produced by OZI forbidden mechanism 17 . 
This leads to a strong suppression (~ 10 -2 ) for heavy quarkonia production 
in hadronic collisions compared to open charm (beauty) production. To pro- 
duce J/ip in the chains by OZI allowed mechanism it is necessary to produce 
2 cc pairs close in rapidity (Fig2.c). Though this mechanism is suppressed 
due to production of extra heavy quark pair it can compete at very high en- 
ergy with the mechanism of single cc pair production. Its contribution can 
be estimated from charm quark fragmentation into heavy quarkonia in e + e~- 
annihilation. Consider production of a cc-pair in q — qq string of Fig. 2. In 
each of q — c and c — qq substrings an extra cc pair can be produced and 
fragment to a given quarkonium state. So it is possible to use an estimate of 
the fragmentation function of c(c) quarks given above or directly experimental 
data from e + e~ to determine a contribution of the corresponding diagrams to 
quarkonia production. This calculation is rather straightforward except of a 
threshold suppression factor. It is clear that at energies of fixed target experi- 
ments y/s = 10 -j- 40 GeV there is a strong suppression for production of J/ip 
and extra DD pair. I shall estimate this suppression factor for an energy of 
HERA-B experiment 18 i?Lah = 920 GeV . Let us denote an extra suppression 
factor compared to suppression of a single cc pair by 7 PP . For its estimation it 
is possible to introduce the same kinematical factors as in e + e _ collisions for 
each subchains q — c and c — qq. For J/ip production at rapidity y=0 7 PP « 0.5. 
Another estimate can be done by assuming that J / ip>- DD system is produced 
by a gluon fusion. This gives j pp ~ 0.4. Using these estimates and taking into 
account that cr^/a 1 ^ w 10~ 2 we obtain that associated production of J/ip with 



charmed hadrons constitute at this energy ~ 10%. At Tevatron energies the 
role of this mechanism is more important and it can (at least partly) explain 
an excess of J / ip production at Tevatron compared to color singlet model. 

For ip' associated production with cc in e+e~ annihilation is not known 
experimentally yet. However its total inclusive yield is close to the one for 
J/ip 19 . If a probabilty of ip' production by c-quark fragmentation is the same 
as for J/ip it will have even stronger impact on ip' production in hadronic 
collision because experimentally for ip' cross section is smaller than for J/ip: 
<J$p/<Jpp ~ 1.6 10~ 3 and associated production can constitute a large fraction 
of the ip' production. 

In conclusion it was demonstrated that the nonperturbativc QGSM model 
predicts a sizable J/ip cc- production in e + e~ annihilation at high energies 
consistent with recent experimental result 6 . In the approach based on 1 /TV- 
expansion in QCD it was shown that a large fraction of cc-quarkonia production 
in hadronic collisions at very high energies can be due to associated production 
with charmed hadrons 
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